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ABSTRACT: We construct anisotropic black brane solutions and analyse the behaviour of 
some of their metric perturbations. These solutions correspond to held theory duals in 
which rotational symmetry is broken due an externally applied, spatially constant, force. 
We find, in several examples, that when the anisotropy is sufficiently big compared to the 
temperature, some components of the viscosity tensor can become very small in units of the 
entropy density, parametrically violating the KSS bound. We obtain an expression relating 
these components of the viscosity, in units of the entropy density, to a ratio of metric 
components at the horizon of the black brane. This relation is generally valid, as long 
as the forcing function is translationally invariant, and it directly connects the parametric 
violation of the bound to the anisotropy in the metric at the horizon. Our results suggest the 
possibility that such small components of the viscosity tensor might also arise in anisotropic 
strongly coupled fluids found in nature. 
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Contents 

1 Introduction 

The AdS/CFT correspondence has emerged as an important tool in the analysis of strongly 
coupled systems, especially for the study of transport properties of such systems. Neither 
analytical nor numerical methods are convenient for calculating these properties on the 
field theory side since they require an understanding of the real time response at finite 
temperature. In contrast, they can be calculated with relative ease on the gravity side, 
often by solving simple linear equations. An important insight which has come out of these 
studies pertains to the behaviour of the viscosity. It was found in KSS [1—3], that for 
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systems having a gravity description that can be well approximated by classical Einstein 
gravity, the ratio of the shear viscosity, 77 , to the entropy density, s, takes the universal 
value 


( 1 . 1 ) 


77 1 

S 4-7T 

This is a small value, compared to weak coupling where the ratio diverges. It was also 
initially suggested that this value is a bound, and the ratio can never become smaller. We 
now know that this is not true [4-7], see also [ 8 , 9], but in all controlled counter-examples the 
bound is violated at best by a numerical factor, and not in a parametric manner. Attempts 
to produce bigger violations lead to physically unacceptable situations, e.g., to causality 
violations, for example, see [10, 11]. However, there is some discussion of a violation of 
the bound in metastable states, see [12]. Also, see [13] for a discussion of violations in a 
superfluid phase described by higher derivative gravity. 

The behaviour of the viscosity discussed above refers to isotropic and homogeneous 
phases, which on the gravity side at finite temperature are described by the Schwarzschild 
black brane geometry. More recently, gravitational backgrounds which correspond to anisotropic 
phases in field theory have also been studied in [ 8 , 14-21] and the behaviour of the vis¬ 
cosity in some of these anisotropic phases has also been analysed, see [22, 23] and [24-29]. 

The viscosity in the anisotropic case is a tensor, which in the most general case, with no 
rotational invariance, has 21 independent components (when the field theory lives in 3 + 1 
dimensions). In [22, 23, 27], where some simple cases were considered, it was found that 
some components of the viscosity tensor can become much smaller, parametrically violating 
the bound in eq.(l.l). For example, in [27], a gravitational solution was considered where 
the rotational invariance of the three space dimensions in which the held theory lives was 
broken from 50(3) to 50(2), due to a linearly varying dilaton . In the solution, the dilaton 
varies along the 2 direction and rotational invariance in the remaining x, y, spatial direc¬ 
tions was left unbroken. The component of the viscosity, called in [27], which measures 
the shear force in the x — y plane, was still found to satisfy the relation, eq.(l.l). However, 
other components of the viscosity did not satisfy it. In particular, it was found that a 
component called 77 _l, which measures the shear force in the x — z or y — z plane, could 
become much smaller, going like 


r]± 87 T T 2 
s 3 p 2 ' 

where T is the temperature and p is the anisotropy parameter. The result, eq.(1.2) is valid 
in the extremely anisotropic limit, when T <C p. A detailed study was also carried out in 
[27] of this extreme anisotropic regime and no instabilities were found to be present. 

In this paper we study many other examples where anisotropic phases arise and show 
that in all of them components of the viscosity can become parametrically small, in units 
of the entropy density, when the anisotropy becomes sufficiently large compared to the 
temperature. Depending on the example, the factor of T 2 in eq.(1.2) can be replaced by 
some other positive power of T. 

A common feature of all our examples is that the breaking of anisotropy is due to an 
externally applied force which is translationally invariant. For example, the linearly varying 
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dilaton considered in [27], and also in section (2) gives rise to such a spatially constant 
forcing function. This follows from the fact that the boundary theory stress tensor is no 
longer conserved in the presence of the dilaton and instead satisfies the equation 

d tl < T> JV >=< 6 > (T<p, (1.3) 

where O is the operator dual to the dilaton, see eq.(6.9) of [27]. Similarly, we consider 
linearly varying axions in section 4.3 and 4.4, and a constant magnetic field in section 4.2. 

Another common feature in our examples is that some residual Lorentz symmetry 
survives, at zero temperature, after incorporating the breaking of rotational invariance. 
Fluid mechanics then corresponds to the dynamics of the goldstone modes associated with 
the boost symmetries of this Lorentz group which are broken at finite temperature. 

In the second half of this paper we give an argument, based on a Kaluza Klein decom¬ 
position of modes, which shows quite generally that in all situations sharing these features, 
in particular where the forcing function does not break translational invariance, appropri¬ 
ate components of the viscosity tensor become parametrically small. These components 
correspond to perturbations of the metric which carry spin 1 with respect to the surviving 
Lorentz symmetry. Let z be a spatial direction in the boundary theory along which there 
is anisotropy and x be a spatial direction along which the boost symmetry is left unbroken, 
then we show that the viscosity component rj xz , which couples to the h xz component of the 
metric perturbation, satisfies the relation, 

VxZ 1 9xx 

s 47T g zz 

where g X x\u=u h , 9zz\u=u h refer to the components of the background metric at the horizon. 
Eq.(1.4) is one of the main results of the paper. It also agrees with the behaviour seen 
in all the explicit examples we consider. This result was first derived for an anisotropic 
axion-dilaton-gravity system in [22], 

In the isotropic case the ratio ^ is unity and we see that the KSS result is 

obtained. However, in anisotropic cases this ratio can become very different from unity and 
in fact much smaller, leading to the parametric violation of the bound, eq.(l.l). 

Note that the result, eq.(1.4), is true for conformally invariant systems, as well as 
systems with a mass gap, when subjected to a constant driving force. Examples of massive 
systems include, for example, gravitational duals of confining gauge theories, [30] and [31]. 
For these cases the temperature should be bigger than the confining scale so that the gravity 
dual is described by a black brane. Also, for some components of the viscosity to become 
significantly smaller than the bound, the anisotropy must be bigger than the temperature. 

Physically a component like r\ xz measures the resistance to shear. For example, if the 
fluid is enclosed between two parallel plates which are separated along the z direction and 
moving with a relative velocity v x along the x direction in a non-relativistic fashion, they 
will experience a friction force due to the fluid, proportional to J) xz d z v x . See Fig (1) and 
the more extensive discussion in section 6 of [27]. Thus the parametrically small values 
obtained here correspond to a very small resistance to shear in anisotropic systems. 


, (1-4) 

U= 
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Figure 1. Picture showing flow of fluid enclosed between two parallel plates separated along the 
^-direction. 


Our results which are quite general, open up the exciting possibility that in nature 
too, strongly coupled anisotropic systems may have a very small value for components of 
the viscosity. It would be very exciting if this behaviour can be probed in experimental 
situations, realised perhaps in cold atom systems, or in the context of QCD. 

This paper is structured as follows. In section 2 we review the earlier discussion of a 
system with one linearly varying dilaton. Some general aspects involved in the calculation 
of viscosity are discussed in section 3. Several examples of anisotropic systems realised in 
gravity are then discussed, including the case with two dilatons in section 4.1, a magnetic 
held in section 4.2, and axions and dilatons, section 4.3 and section 4.4. The general 
argument based on a Kaluza Klein truncation is given in section 5. We end with conclusions 
in section 6 . The appendices A, B and C contain additional important details. 

2 Brief Review of The System With One Dilaton 

Here we briefly summarise some of the key results in [27] which considered a linearly varying 
dilaton cf> = p z in asymptotically AdS 5 , for a theory with action 
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Here 2 k 2 = 167 tG is the gravitational coupling and G is the Newton’s Constant in 5- 
dimensions. At zero temperature the near horizon solution was found to be AdS 4 X R, 


ds 2 = ~^u 2 dt 2 + + \ u 2 (dx 2 + dy 2 ) + dz 2 . ( 2 . 2 ) 

o -^u z o o 

The radius of AdS 4 , R\ = 3/4, in units where A = 1 . We see in eq.( 2 . 2 ) that the metric 
component g zz becomes constant due to the extra stress energy provided by the linearly 
varying dilaton. The AdS4 x R solution is in fact an exact solution to the equations of 
motion. 

At small temperature, T <C p, the geometry is that of a Schwarzschild black brane in 
AdS4 x R. The viscosity is related, using linear response, to the retarded two point function 
of components of the stress tensor, and the latter, using Ads/CFT, can be calculated from 
the behaviour of appropriate metric perturbations in the bulk. The answer for rj xy , which 
is denoted as rj 1 | and for r] xz ,r] yz , which are equal and denoted as r/j_, is given in eq.(2.3) 
and eq.(2.4) below: 


m = _l 

S 4-7T ’ 

rj± 8 irT 2 

s 3/o 2 


(2.3) 

(2.4) 


with s being the entropy density. 

We see that g± in units of the entropy density becomes parametrically small in the limit 
of high anisotropy. The fluid mechanics in this high anisotropy limit was also systematically 
set up in [27] and it was shown that, as expected, this small viscosity component results in 
a very small shear force on two suitably oriented parallel plates which are moving with a 
relative velocity and enclose the fluid. 


3 More Details On The Calculation Of Viscosity 

Before proceeding it is worth giving some more details on the calculation of the viscosity 
for the one dilaton system above. These features, as we will see, will be shared by all the 
examples we consider subsequently in this paper. The analysis that follows will also reveal 
the central reason for why the viscosity in units of the entropy density can become so small 
in anisotropic systems. 

With anisotropy, the viscosity is a tensor, rjijki, in general with 21 components. Using 
the Kubo formula these can be related to the two point function of the stress energy tensor 
as follows , 

Vij,kl = - lim - Im[Gj H (w)], (3.1) 

(jJ —^0 UJ J 

where 

Gg k i(u>,0) = J dt dxe^ 6(t) ([Tijit, x),T w ( 0 , 0 )]), (3.2) 

and Im denotes the imaginary part of the retarded Green’s function. 

From the AdS/CFT correspondence the two point function of can be calculated in terms 
of the behaviour of metric perturbations, and in this way the viscosity can be obtained. 
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In the one dilaton system considered in section 2 , the solution has an SO( 2) rotational 
invariance in the x — y plane, as is evident from the metric (2.2). For simplicity we denote 
the r) X z, xz component as r) xz , and g yz , yz as r] yz etc. Due to the £0(2) invariance we get that 
Vxz = Vyz = V±" These components are related to the behaviour of the h xz , h yz components 
of metric perturbations, which carry spin 1 with respect to £0(2) symmetry. 

We now proceed to introduce the h xz perturbation in the metric as follows 


ds 2 = - gtt{u)dt 2 + g uu (u)du 2 +g xx (u)dx 2 + g yy dy 2 + g zz {u)dz 2 

+ 2 e~ Iujt Z(u)g xx (u)dx dz, 


(3.3) 


where Z(u) is the required perturbation of interest. We can show that the other modes 
decouple from Z(u ) and hence we can consistently set them to zero. Here we follow closely 

[32]- 

One finds that the mode Z(u ) obeys an equation of the form 

du (V~ 9 P(u)g uu d u Z (u)) - oj 2 N(u)g tt Z(u ) = 0, (3.4) 

The functions P(u), N(u) are given in terms of the background metric , with 

P(u) = g zz g xx ■ (3.5) 

In effect, eq.(3.4) arises from an action 

S=-j - lN(u)g“(a,Z) 2 ] (3.6) 

(we are neglecting the dependence on the spatial x l coordinates here). Using AdS/CFT we 
can find the response in terms of the canonical momentum 


n(rt, lu) = — 


1 


167tG 


V^P(n)g uu d u Z(u). 


(3.7) 


The retarded Green’s function is then given by the ratio of the response over the source, 

n(«, u 


G ret = - 


Z(u, uj) 


(3.8) 


leading to the result from eq.(3.1) 


f]± = lim T 


n(u, u) 


o iuZ(u, oj) 


(3.9) 


We now show that the RHS of eq.(3.9) can also be evaluated near the horizon, u = uh, 
instead of u —> oo. Since we are interested in the limit w-^Owe can neglect the second 
term in eq.(3.4) leading to 

<9 u n = 0 (3.10) 

upto 0(uj) 2 . This gives 

n = C, (3.11) 








where C is independent of u . Next, it is easy to see that there is a solution of eq.(3.4) 
in the u —» 0 limit in which Z is simply a constant. This solution also meets the correct 
boundary condition at u —> oo, since, as can be seen from eq.(3.3), the non-normalisable 
mode must go to a constant at u —> oo. Putting all this together we find that to leading 
order in the u —> 0 limit both II and Z are constant and thus the ratio in eq.(3.9) being 
independent of u can also be evaluated at the horizon. 

As a result we get 


rj _|_ = lim 


II(«, oj) 


oj->o iuj Z(u, oj) 




(3.12) 


Demanding regularity at the future horizon , we can approximate the behaviour of Z 
as follows 

Z ~ e ~ ia, ( t+r *), (3.13) 


where r* is the tortoise coordinate, 


r* 



(3.14) 


It then follows that 


The entropy density is 


1 


V± = 


16nG 


P(u H ) 


-g 


QttQuu 


U—^Uh 


1 yf^g 
4G yjg uu gtt 


(3.15) 


(3.16) 


to 


Using the value of P(u) from (3.5) and using eq.(3.15) and eq.(3.16) this finally leads 


?7-L 

s 


i gxx 
4vr g zz 


U H 


(3.17) 


We now see why anisotropic systems will generically be different from isotropic ones. 
For an isotropic system rotational invariance makes the ratio ^ = 1, leading to the KSS 
bound, eq.(l.l). However in the anisotropic case in general this ratio will not be unity and 
thus the ratio of 77 /s can become smaller than In the one dilaton system this is what 
happens leading to the result, eq.(3.17). In the rest of this paper we will find many more 
examples of this type, where anisotropy will allow different metric components to shrink 
at different rates and attain different values at the horizon, thereby leading to violations of 
the KSS bound. 


4 Additional examples with anisotropy 

4.1 Anisotropic solution in two dilaton gravity system 

To generalise the example in section 2, we consider next the case of gravity, with a negative 
cosmological constant, two massless scalar fields, <f>i and (j )2 ; both of which we now call 
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dilatons, in 5 spacetime dimensions with action, 






(4.1) 


Both the dilatons are turned on to be linearly varying, but along different directions: 


(4.2) 


01 = piy, 02 = P2Z. 


The zero temperature near horizon solution is now given by AdS 3 x R . x R (we have set 


A=l): 



(4.3) 


We see that there are now two different mass scales, pi, P 2 which characterise the 
anisotropy. In appendix A we show that this near horizon geometry interpolates smoothly 
to asymptotically AdS§. The 50(2,2) symmetry of AdS 3 is preserved all along this inter¬ 
polation. 

At small temperature, T <C pi, P 2 , the near-horizon solution is given by : 



The horizon lies at 


(4.5) 


u = Uh = 7 tT. 


The computation of the shear viscosity follows the discussion in [27] quite closely. The 
near-horizon AdS 3 has 50(1,1) Lorentz invariance in the t, x directions. The metric pertur¬ 


bations can be classified in terms of different spins with respect to this 50(1,1) symmetry. 
The viscosity component r] xz , given by, 



(4.6) 


can be calculated by considering a metric perturbation Z(u) defined so that the full metric 
with the perturbation takes the form, 


ds — gtt(u)dt T g uu {u)du T g X x( , u)dx g-yydy“ gzzid^)dz 

_ 1 _ 2 e~ lwt Z(u)g xx (u)dxdz. 


(4.7) 


This component has spin 1 with respect to the 50(1,1) symmetry. It turns out that result¬ 
ing analysis is quite similar to that in section 3 and this perturbation satisfies an equation 
of the type given in eq.(3.4), with P(u) given by eq.(3.5). The conjugate momentum II is 
also given by eq.(3.7) with P(u) given by eq.(3.5). As a result y xz is given by eq.(3.15). 

The entropy density is given by 


1 V~9 


40 yjg U u9tt 
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s = 


(4.8) 








This gives, 


which using eq.(4.4) becomes 


Similarly, for r] xy we get 


Pxz 1 9xx 

S 47 t g zz 

r]xz _ 47 tT 2 

s pI 


Vxy 1 9xx _ 4ttT 

S 47r g yy p\ 

U H 


(4.9) 


(4.10) 


(4.11) 


We see from eq.(4.9), eq.(4.11) that the relative ratio of y/s for these components is 
determined by the ratio of the metric components as one approaches the horizon. 


4.2 Viscosity in the Presence of a Uniform Magnetic Field 

Here, for completeness, we briefly review a situation where the anisotropy is generated due 
to a magnetic held which has been studied in considerable depth in [28]. We refer to [28] 
for details. We start with a system with the action 

S = J d 5 x^g{R +12 A- ^F 2 ), (4.12) 

and consider a solution where the magnetic held 


F — B 

1 yz — 1J ) 


(4.13) 


with B being a constant. Such a system was also considered in [33]. 

The resulting near horizon solution at zero temperature is now again AdS 3 X R X R, just 
as in the two dilaton system, with rotational invariance also preserved in the yz plane. 
The metric is (we have set A—1) 

ds 2 = —3 u 2 dt 2 H— \^du 2 + 3 u 2 dx 2 H- ^=\B\dy 2 H- \B\dz 2 . (4.14) 

3 ir 2y3 2y3 


The radius of AdS%, R 2 = 1/3, in units where A = 1. 

At small temperature, T <C B the solution is a black brane in AdS 3 x R X R with 
metric 


ds 2 = — 3w 2 (1- y)dt 2 + 


1 


1 


r 


3« 2 (1-A,) 


dir + 3 irdx 2 H- ^=\B\dy 2 H- -^\B\dz 2 , (4.15) 


2V3 


2 ^ 


where c is given in terms of T as follows 


c = 


4tt 2 T 2 

' 


U = u h = ^ttT. 


(4.16) 
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The horizon lies at 


(4.17) 












The viscosity components r] xy = 7 ) xz = rj±. To calculate rj± we consider the h xz 
component of metric perturbation, so that the full metric is of the form 

ds 2 = -g tt {u)dt 2 + g rr (u)dr 2 +g xx (u)dx 2 + g yy {u)dy 2 + g zz {u)dz 2 

. (4.18) 

+ 2e lui Z(u)g xx (u)dxdz, 

with Z(u ) being the perturbation that we need to study. One can easily show that the 
other modes decouples from Z{u ) and so can be consistently set to zero. 

We find that the resulting analysis is again quite similar to that in section 3 . This 
perturbation satisfies an equation of the type given in eq.(3.4), with P(u ) given by eq.(3.5). 
The conjugate momentum II is also given by eq.(3.7) with P(u ) given by eq.(3.5). 

The resulting value for the viscosity is given by 


V± 

s 


1 9xx 


4vr g zz 


U H 


Substituting the metric components from (4.15) above we get that 


(4.19) 


V± = 

s V 3 n \B\ 


(4.20) 


As discussed in [28], this example may be relevant in the study of QCD, perhaps for 
heavy ion collisions, and also in the core of neutron stars where strong magnetic fields can 
arise. 


4.3 The Dilaton-Axion System 

In the examples considered so far, the near horizon geometry was of the form, AdS X R n , 
with the metric components along the R n directions not contracting as one gets to the 
horizon. It is worth considering other situations where the near horizon geometry is of 
Lifshitz type instead, with metric components along all the directions contracting as one 
approaches the horizon but at different rates. 

An easy way to construct such an example involves a system consisting of gravity with 
an axion and dilaton with action, 

Sbuik = J d b xy/^g (r + 12A - ^( dc,b) 2 - (<9y) 2 ) > ( 4 -21) 

containing the parameter a which enters in the dilaton dependence of the axion kinetic 
energy term. Earlier work in [22] considered the case with a = 1. The case a = — 1 has 
SL( 2, R) invariance. 
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It is easy to see that by turning on a linear profile for the axion one obtains an extremal 
solution whose near horizon limit is given by ( setting A—1) 


du 2 


ds = R —u~dt H-y- + u~dx + udy + p u 1 + 2q2 dz/ 


u* 


ci = 


X = ci p z, 

^2(3 + 8a 2 ) 
(l + 2a 2 ) ’ 
2 3 + 8cr J 

R = 4 + 8a 2 ' 


(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 


This solution breaks rotational invariance along the z direction due to the linearly 
varying axion, and p is the mass scale which characterises this breaking of anisotropy. 
We see that all components of the metric along the spatial directions now shrink as one 
approaches the far IR, but the rate at which the g zz component vanishes is different from 
the other spatial components, g xx ,g yy - Let us also note that for a = 1 the solution above 
agrees with [15]. 

At small temperature T -C p the resulting solution has a metric given by 


ds 2 = R 2 [ —u 2 f{u)dt 2 + 


du 2 

u 2 f(u] 


+ W'dx + u dy" + p u 1 + 2a2 dz " 


(4.27) 


where R 2 is as given in eq(4.26) above and f(u) is given as 


1 - 


/l6vrT\ p 
\ p 2 u ) 


(4.28) 


where p = . The axion continues to be linear as in the solution eq.(4.23) and the 

dilaton is given by eq.(4.24). 

The horizon in eq.(4.27) is at 


167tT 


u = u h = 




(4.29) 


Let us now turn to computing the viscosity. The shear viscosity component rj xy satisfies 
the KSS bound in eq.(2.3) . Next consider the component rj xz = r) yz . To compute this 
component we can consider the h xz component of metric perturbation, so that the full 
metric is of the form 


ds~ — gtt{u)dt~ + g uu {u)du ~\~g xx {u)dx T gyy(u)dy T g zz {u)dz 

+ 2 e~ lult Z(u)g xx (u)dxdz, 


(4.30) 


where Z(u ) is the perturbation that we need to study. The dilaton and axion are unchanged 
and are given by eq (4.24) and eq (4.23) respectively. One can easily show that the other 
modes decouples from Z(u ) and so can be consistently set to zero. 
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We again find that resulting analysis is similar to that in section 3 and the perturbation 
satisfies an equation of the type given in eq.(3.4), with P(u) given by eq.(3.5). The conjugate 
momentum II is also given by eq.(3.7) with P(u) given by eq.(3.5). As a result r] xz is given 
by eq.(3.15). 

Thus, substituting the metric components for the finite temperature solution (4.27) we 


get 


771 1 q xx , T . 2 

d± = _^ (_) 1 + 2 Q 2 _ 

S 4-7T g zz p 


(4.31) 


The dependence on T in eq.(4.31) follows from the metric eq.(4.27) and the dependence 
on p is then obtained on dimensional grounds. Let us note that the temperature T which 
appears in eq.(4.28) could be related to the temperature as measured in the asymptotic 
AdS coordinates by a rescaling. By the asymptotic AdS coordinates we mean those in 
which the metric takes the standard form: 


ds 2 


— u 2 dt 2 + + u 2 (dx 2 + dy 2 + dz 2 ) , 

u- 


(4.32) 


This is also true for the x,y coordinates in eq.(4.27) and the corresponding coordinates 
which appear in eq.(4.32). and also for the z coordinate in eq.(4.27) which is related to the 
corresponding coordinate in eq.(4.32) by a p dependent rescaling in general. These rescaling 
factors have to be determined if the coefficient in eq.(4.31) is to be fixed. To do so, one 
needs to find the full interpolating geometry from the near horizon region, described by 
eq.(4.27) , to the asymptotic AdS region, eq.(4.32). 

We have carried out such a numerical interpolation for a = ±1, for which, eq.(4.31) 
becomes, 

— ~ (-) 2/3 . (4.33) 

s p 

We find, within the accuracy of our numerical calculation, that there is no rescaling of the 
T,x,y coordinates while the 2 coordinate is rescaled by a non-trivial p dependent factor. 
One consequence is that the temperature T which appears in eq.(4.31) is the same as the 
temperature measured in the held theory. 


4.4 The two Axion-one Dilaton System 

For good measure, as another example, we consider a system consisting of gravity with two 
axions and one dilaton described by the action 

Sbuik = 7^2 J d b Xy/^g (r + 12 A - ^{d^f - ^e 2 a 4 > (dxi ) 2 ~ ^e 2a<t> ( 8 x 2 )^ ■ (4.34) 

In this case we will see that for a suitable profile for the two axions, the AdS 4 sym¬ 
metry of the near-horizon geometry is broken further to AdS 3 , with now two of the spatial 
directions, y, z. being characterised by non-trivial Lifshitz exponents. 
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The linear profiles for the two axons and resulting near horizon solution is given by 
(setting A—1) 


ds z = R z ( — u 2 dt 2 + + u z dx z + p z «!+ 4q2 dy z + p z u x + 4a ' 2 dz z ) , 


du 2 


OO O S °T, O O 8 Ay o 

I ^ » . t 1 rl .v'7 ,7, , t „ , 1 




c = 


1 + 4a 2 


Xi = c p y, 

X2 = C P Z, 

4 alog(ii) 

1 + 4a 2 ’ 

\/l + 8ct 2 , 


^2 _ 1 + 8a " 


2 + 8a 2 ' 


(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 


This metric in this solution has AdS^ invariance, and also a scaling symmetry under 
which y,z transform with a non-trivial exponent. The linearly varying axions break this 
scaling symmetry, and also the rotational invariance along the y and z directions, with p 
being the mass scale which characterise the breaking. 

At small temperature T <C p the resulting solution has a metric 

du 2 
u 2 f(u ) 

where R 2 is as given in eq.(4.40) above and f(u) is given as 

16 vrTy 
p 2 u J ’ 




R 2 ( — u 2 f(u)dt 2 + 


+ u 2 dx 2 + p 2 u !+4“ 2 dy 2 + p 2 a 1 +-*“ 2 dz 2 ^j , (4-41) 


where p = • 

The two axions continue to be linear as in the solution eq.(4.36), eq.(4.37) and the 
dilaton is given by eq.(4.38). 

The horizon in eq.(4.41) is at 


u = u h 


167 tT 


(4.43) 


The r) xy and q xz components of the viscosity are the same., we denote them by r/y. To 
calculate these components we consider the h xz component of metric perturbation, so that 
the full metric is of the form 


ds — gtt{u)dt~ + g uu (u ) du" +g xx (u)dx T gyy(u)dy g zz (u)dz 

_l_ 2 e~ lujt Z(u)g xx (u)dxdz, 


(4.44) 


where Z{u) is the perturbation that we need to study. 

The dilaton and axions are unchanged and are given by eq (4.38) and eq (4.36) , eq (4.37) 
respectively. One can easily show that the other modes decouples from Z(u) and so can be 
consistently set to zero. 
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As in the previous cases , the analysis here is similar to that in section 3 and this 
perturbation satisfies an equation of the type given in eq.(3.4), with P(u) given by eq.(3.5). 
The conjugate momentum II is also given by eq.(3.7) with P{u) given by eq.(3.5). As a 
result r] xz is given by eq.(3.15). 

Thus, substituting the metric components for the finite temperature solution (4.41) we 


get 


V± _ ^ l+4o 2 

S 47 t g zz \pj 


(4.45) 


For the case a = ±1, eq.(4.45) becomes, 


V± 

s 



(4.46) 


Interestingly, both in eq.(4.31) for the one axion case, and in eq.(4.45) above we see 
that the maximum value the exponent governing the temperature dependence can take is 
2, and the minimum value, for a = oo, is 0. 


5 Kaluza Klein Reduction 


The previous sections dealt with a number of examples where anisotropic situations gave 
rise to small values for the viscosity to entropy ratio. One common feature of all these 
examples was that the breaking of isotropy was due to a spatially constant driving force. 
For example, the dilaton considered in section 2, gives rise to a force proportional to the 
gradient of the dilaton which is a constant since the dilaton varies linearly. One way to see 
this is by noting that the stress tensor is no longer conserved and satisfies the equation 

d tl < V‘ v > = <6> <9>, (5.1) 


as discussed in eq.(6.9) of [27]. Similarly, we consider linearly varying axions in section 4.3 
and 4.4, and a constant magnetic held in section 4.2. 

In this section we will present a general argument which should apply to all such situa¬ 
tions where the breaking of isotropy occurs due to matter fields which give rise to a spatially 
constant driving force. We will also assume that a residual AdS symmetry is preserved in 
the bulk, and a corresponding Lorentz symmetry is left intact in the boundary theory. Fluid 
mechanics then corresponds to the dynamics of the goldstone modes associated with the 
boost symmetries of this Lorentz group. The components of the viscosity which give rise 
to the violation of the KSS bound in the examples considered above correspond to metric 
perturbations which have spin 1 with respect to the surviving Lorentz symmetry. Let z 
be a spatial direction in the boundary theory along which there is anisotropy and x be a 
spatial direction along which the boost symmetry is left unbroken then we will present a 
general argument below showing that the viscosity component r] xz , which couples to the 
h xz component of the metric perturbation, satisfies the relation, 


Vxz 

S 


1 9xx 

47r g zz 


u=u h 


(5.2) 
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where g X x\u=u h , 9zz\u=u h refer to the components of the background metric at the horizon. 
Eq.(5.2) is the main result of this section and one of the main results of this paper. We note 
that it also agrees with all the examples considered above. This result was first obtained 
for an anisotropic axion-dilaton-gravity system in [22], An analysis using RG flow and KK 
reduction, for this system, was carried out in [26] along the lines of [2, 32], 

For a case with a residual AdS^+i factor in the metric, the basic idea behind the 
general analysis will be to consider a dimensionally reduced description, starting from the 
original D + 1 dimensional theory and going down to the AdSd+i space-time. Different 
Kaluza Klein (KK) modes in the extra dimensions will not mix with each other since the 
effects breaking rotational invariance are in effect spatially constant. For example, for cases 
where there are linearly varying fields, like axions or dilatons, this will be true since the 
equations of motion involve only gradients of these fields which are spatially constant. The 
non-mixing of the KK modes will greatly ease in the analysis, since we can use the standard 
formulae of KK reduction and moreover truncate the analysis to the zero modes in the extra 
dimensions. The off diagonal components of the metric, whose perturbations carry spin 1 
and which are related to the viscosity components of interest, will give rise to gauge fields 
in the dimensionally reduced theory. By studying the conductivity of these gauge fields, 
which can be related easily to the spin 1 viscosity components we will derive the result in 
eq.(5.2). 

The study of more complicated situations where the breaking of rotational invariance 
is due to a driving force that also breaks translational invariance is left for the future. 


5.1 The Dimensionally Reduced Theory 

To start, we will consider the case where D = 4 and d = 3, so that a residual AdS 4 
symmetry survives, and the asymptotic geometry, towards the boundary, is AdS$. In this 
case we start with 5 dimensions with a gravitational action : 


S = ^ J d 5 X ^j (R+ 12 A). (5.3) 

Here 2 k 2 = 167 tG is the gravitational coupling with G being Newton’s Constant in 5- 
dimension and we set A=l. 

Parametrising the 5 dimensional metric by 


( er^ u) g^ + e 2 ^A^ \ 

^ e 2 ^A u e 2 ^ ) ’ 

and taking all components to be independent of the z direction which we take to be the 
compactification direction, gives 

S = 2k? / “ TfW) 2 - Tf f2 + , (5.5) 

where we have dropped total derivatives . 

We also note that in our choice of parametrisation , 



9zz = 




(5.6) 
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The coefficient of the first term in the matrix in eq.(5.4) was taken to be e - ^ so that the 
resulting 4 dimensional action is in the Einstein frame, k which appears above is related 
to the 5 dimensional gravitational coupling k by 

L _ 1 
2 k? ~ 2^’ 

where L is the length of the compactified 2 direction. 

So far we have neglected any matter fields. Consider for concreteness the case of the 
axion-dilaton system considered in section 4.3 with action eq.(4.21) with a— 1. Inserting 
the background solution for the axion 


(5.7) 


X = a z, (5.8) 

and taking the dilaton to be independent of z we get from the kinetic energies of the dilaton 
and axion, 


s=T / „v=5 (-^ - \m 2 - ■ (*■») 

We see that there is an extra term which depends on the gauge field and which gives rise to 
a mass for it. This term arises due to the linearly varying axion, eq.(5.8) and is tied to the 
breaking of translational invariance due to this linear variation. We see that the terms in 
eq.(5.5) and eq.(5.9) involving the gauge field are quadratic in this field and can be written 
as 

s = i?j dlx ^ s {jisu) p2 - i™ 2 *")- 42 ) • < 5 - 10 > 

where 

m 2 (u) = 2 a 2 e 2 ^ u \ (5-11) 

and 

2 1 ( ' = e 3 ^ = {g zz (u))% . (5.12) 

9es( U ) 

The solution in the near horizon region for this dilaton-axion system was given in 
eq.(4.22) with a— 1. It is easy to see from this solution that 

= p 3 u 2 , (5.13) 

9eS 

and 44 

m 2 (u) = —p 2 uz, (5-14) 

and therefore that the gauge coupling and mass vary with the radial coordinate. 

Similarly, in other cases where there is also a breaking of translational invariance we 
will get both a kinetic energy term and a mass term, and in general both the gauge coupling 
and the mass will vary in the radial direction. For the subsequent analysis we will analyse 
the perturbations of the gauge fields in the 4 dimensional theory given in eq.(5.10). Such a 
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system was considered in [34, 35] and our subsequent discussion closely follows this reference. 
As we will see later, the conductivity of these gauge fields can be related easily to the spin 
1 viscosity components using which we will derive the result in eq.(5.2). Let us mention 
for now that the essential reason for this is that the two-point correlator of the current 
operator gives the conductivity of the gauge held, while the two-point stress tensor in the 
higher dimensional theory is related to the viscosity. Since the gauge held is obtained from 
the spin 1 component of the metric in the higher dimensional theory, these two correlators 
are closely related. 

The 3 + 1 dimensions, include time, t, the radial direction u, and additional space direc¬ 
tions, one of which we denote by x. To study the conductivity we consider a perturbation 
for the x component of the gauge held, 


A x (x,t,u ) 


dud 3 k 

(2vr) 4 


e -^t+k.x Z ( U ,u). 


(5.15) 


This gauge held perturbation decouples from the rest (we have set perturbations of the 
axion to vanish even before the KK reduction in the example above, this turns out to be a 
consistent thing to do). Z(u, oj) satishes the equation 


d_ 

du 


{ N (u) — Z ( u ,u)) 


lo 2 N{u) g uu g U Z(u, ur) + M(u)Z(u, u>) = 0, 


(5.16) 


with 


and 


N{u) = V^^g xx g uu , 
9eS 


M(u) = - 


m 2 (u)y/=g 
2 gxx 


(5.17) 

(5.18) 


Treating the radial coordinate u as the analogue of time we can read off the “momentum” 
conjugate to Z from eq.(5.10) to be 


II(u, u ;) 


5S 

5Z'(u, —uj) 


1 

2 ^ 


N(u)Z'(u,uj), 


(5.19) 


where Z' = 4^Z(u,uj) and N(u) as given in eq(5.17). 
The conductivity is given by 


cr(u , ui) 


n(rt, w) | 

iuZ(u,uy v ^ 00 ’»->o’ 


(5.20) 


where Z and II are the asymptotic values of the perturbation and conjugate momentum 
defined in eq.(5.19) in the region u -+ oo. 

We assume that the underlying higher dimensional geometry is asymptotically AdS§ 
space and that the back reaction due to the matter fields which break the rotational in¬ 
variance dies out compared to the cosmological constant in this asymptotic region. This is 
true in all the examples studied above where the geometry becomes AdS§ when u —> oo. 
It is then easy to check, as discussed in appendix B that the ratio on the RHS in eq.(5.20) 
becomes independent of u when u —> oo. 


17- 






We can write a(u, uj) as the sum of real and imaginary parts as 5?(cr(«, ui))+i Q(a(u, a;)). 
We will be interested in the real part 3?(<r) since that is related to the viscosity components 
of interest. It is easy to see from our definition, eq.(5.20) that 


3?(cr(it, w)) 


/ II(n, oj)Z(u, — uj) V 
\^ u jZ(u,uj)Z(u,—uj) 


(5.21) 


where n(ti, w) is defined in eq.(5.19). 

To evaluate the RHS in the limit cu —>• 0, it will be sufficient to consider the leading order 
behaviour of the denominator. Since Z(u, uj) is real to leading order when w->0we obtain 


5?(ct) 


9(II(u, uj)Z(u, —uj)) | 

UJ z 2 (u) w - s -0' 


(5.22) 


The numerator of RHS of eq.(5.22) is independent of u (appendix C) and can therefore 
be evaluated at u = Uh instead of u —> oo. After some more simplification this gives 


5ft(<j) = (7 H 


( Z(u h ) \ , 

yz(u oo) J 


(5.23) 


where an is the conductivity evaluated at the horizon and its expression is given by, 




1 


2 K 2 g 2 eS 


U=Uh 


(5.24) 


See appendix C for more details. 

To proceed we need to evaluate the ratio z ^u^oo) • ^ or P ur P ose we g° back to 
the underlying higher dimensional theory with which we started in which the gauge field is 
actually an off diagonal component of the metric, eq.(5.4). The background about which 
we are calculating the behaviour of the perturbation is diagonal in the metric with all 
components being only a function of u. Now consider a coordinate transformation x —> 
x + az, with all the other coordinates remaining the same. It is easy to see that under this 
transformation the metric now acquires an off-diagonal component 


$9xz — O (jxx ■ 


(5.25) 


with all the other components of the background metric staying the same. Note that in our 
notation the hatted metric refers to the 5 dimensional one while the unhatted metric refers 
to the 4 dimensional Einstein frame metric, see eq.(5.4). 

Since we have merely carried out a coordinate transformation it is clear that 5g xz in 
eq.(5.25) must satisfy the equations of motion for small perturbations about the starting 
background. Comparing with eq.(5.4) we find that this corresponds to turning on a gauge 


held 


A~ = a 


9xx 


which must therefore solve the equation (5.16) in the limit uj —» 0 with 


(5.26) 








In this way we can exploit the co-ordinate invariance of the underlying higher dimensional 
theory to obtain a solution for Z(u) in the co —> 0 limit. More over it is easy to see that this 
solution meets the correct boundary condition at u —>■ oo. As was mentioned above, we are 
assuming that the higher dimensional metric is asymptotically AdS§ space. The ratio 
therefore goes to unity and Z{u ) goes to a constant which is the correct behaviour needed, 
as is also discussed in appendix B. 

With the solution eq.(5.27) at hand we can now evaluate the ratio z^v^lo) • The 
arbitrary constant a drops out and we get that 

Z(uh) _ (jxx 

Z(u —>• oo) e 2 ^ 

Substituting in eq.(5.23) and using eq.(5.24) we get that the conductivity is given in 
terms of g 2 ff and various metric opponents at the horizon by 

1 f 9XX \ 

ff,Is w*) 

From eq.(5.12) and using eq.(5.6) from our parametrisation eq.(5.4) , we finally get that 


(5.29) 

11 = 


(5.28) 

U= 


a = 


1 Q 2 
vxx 

2k 2 y/gTz' 


(5.30) 


Note that we have been able to obtain an expression independent of m 2 that only 
depends on the metric components g X x,9zz hi the 5 dimensional theory. In the subsequent 
discussion we somewhat loosely denote Re(a) by a itself. 


5.2 The Viscosity To Entropy Ratio 

The next step is to relate the conductivity obtained above to the viscosity. This is in fact 
straightforward. Kubo’s formula relates the components of the viscosity to the two point 
function of corresponding components of the stress tensor in eq.(3.1). This two point 
function is obtained by calculating the response to turning on suitable metric perturba¬ 
tions in the bulk. We will be assuming, as was mentioned above, that asymptotically the 
background metric is AdS§. Thus as u —> oo, g^ u —> u 2 5^ l , for all components other than 
along the u direction, as discussed in appendix B. The off - diagonal metric perturbations 
required for the shear viscosity then behave like 

^9fj,u = n 

as u —> oo, where is independent of u. The viscosity component rj xz is then given by 

dxz = ~ T xz (ki,w)T xz (k2,uj) > ^ (5.31) 

where the prime subscript on the RHS means that the overall energy momentum conserving 
delta function has been removed. From AdS/CFT we have that 

- -> S 2 S 

< T xz (ki)T xz (k 2 ) > = - —=r~ -=r-. (5.32) 

5hx Z (h)5h xz (k 2 ) 
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The conductivity in an analogous way is given by 


a = -—3 (< J x (ki,u)J x (k 2 ,u) >') 

(jJ \ / 


fci } fc2—>0,o;—>0’ 


which in turn can be calculated from the bulk response since 

< J x (h)J x (k 2 ) > = 


5A x (ki) SA x (k 2 ) 


(5.33) 


(5.34) 


On comparing with eq.(5.4) we see that the zero mode of h zx in the z direction is in 
fact A x . This shows that r/ xz and a are essentially the same upto one minor factor of L 
the size of the 2 direction. This factor arises because the prime subscript in eq. (5.31) and 
eq.(5.33) are different, in the first case the momentum conservation delta function removed 
includes a delta function in the z direction, whereas in the case of the conductivity it does 
not include this delta function. Accounting for the difference gives 

<j 

Vxz = (5.35) 


The entropy density in the 5 dimensional theory is given by 


2 ir 2 n p —-—-— 
£ = = I? y/ 9xx9yy9zz > 


K 


K 


(5.36) 


(this is also the same as the entropy density in the 4 dimensional theory divided by L). 
From eq.(5.35) , eq.(5.29), eq.(5.6), eq.(5.36)and eq.(5.7), we can now write the ratio 


Vxz _ 
s s 


1 


9eS 



^ \J9xx9yy9zz 


u=u h 


(5.37) 


Using eq.(5.12), eq.(5.6) in the above expression and using isotropy along x and y, we arrive 
at the following result 

VXZ _ 1 9XX I [r 

s 47r g zz ' u=u h 

This general result agrees with the ones we obtained in all the examples we studied in the 
previous sections. We see that independent of the details of the matter fields which were 
responsible for the breaking of the rotational symmetry we get a general result in eq.(5.37). 
This result shows that when the ratio of the metric components ^ at the horizon becomes 
smaller than unity the KSS bound will be violated. 


5.3 Generalisation To Case with Additional Directions 

In the preceding discussion of this section we have considered the dimensional reduction 
from 5 to 4 dimensions. However, it is easy to generalise these results for the case where we 
start with D +1 dimensions and KK reduce to d +1 dimensions. In fact, this generalisation is 
needed for the situation discussed earlier with a magnetic field where the residual symmetry 
arises due to an AdS 3 factor instead of an AdS 4 in the geometry. Our analysis closely follows 


20 







[36]. The dimensional reduction in this case will give rise to D — d gauge fields . 
Following [36], we parametrize the higher dimensional metric as : 


Fb A 1 ) 


( 9ab ) = 


g^ + A^A 
4 (1) 

■**-VOL 


I /7 



where the D + 1 dimensional vielbein is given by 


e r , = 


= v 



(5.39) 


(5.40) 


with G a p = E^dabEp and = e()r/ rs e*. Here a,/3 = 1, ..D — d denote the directions over 
which the reduction has been carried out and //, v = 0,1, • ■ ■ d are the ones left in the lower 
dimensional theory. It also follows from the parametrisation that 


\f~9 = 


(5.41) 


where G is the determinant of the internal metric G a g. Additional matter fields required 
for breaking rotational invariance which also break the translational invariance in the com- 
pactihed directions give mass terms for the gauge fields, which will vary in general in the 
radial direction. Neglecting these additional matter fields for now we start with the action 

fd D+1 x^{R + A] 

As shown in [36] the dimensionally reduced action in d + 1 dimensions becomes 
d d+l xy/^g (^R + A + (f^d^dyd + ^g^d^G a/3 d u G al3 
-\g» p g vX G af) F$ a F^ , 

where 

cj) = ~^log det {G a p) => e -0 = VG, 
where G is the determinant of the internal metric G a p, 



(5.42) 


(5.43) 


F^ a = d^ a -d u A^ a , (5.44) 

and k which appears above is related to the 5 dimensional gravitational coupling k by 

L D ~ d _ 1 
2 k 2 ~ 2 ^’ 

where L D ~ d is the volume of the compactihed directions . 

For simplicity we assume that the internal metric G a p is diagonal and focus on the g xz 
component of the metric perturbation (where x represents a spatial direction along which the 
boost symmetry is left unbroken and z represents an anisotropy direction in the boundary 
field theory). Comparing the last term in the action eq.(5.42) with the kinetic energy term, 
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y/—g ^ ^ -F 2 ^, as given in eq.(5.10), we then find the effective gauge coupling, for the 

corresponding gauge field A is 

= e~^g zz - (5.46) 

9 is 

As mentioned above, additional matter fields give rise to mass terms for the gauge 
fields. We will also take these mass terms to be diagonal for simplicity. The resulting 
equation for the x component of the gauge field A x is then of the form given in eq.(5.16), 
where we have expanded A x as given in eq.(5.15). It can then be argued (see Appendix C 
for details) that the conductivity in the lower d+ 1 dimensional theory 1 is given by 


= ^2 I 


1 

2 k? 


Quu 

9tt 


N(u ) 


u=u h 


Z(u h ) 
Z(u s oo) 


guii , - l X x v 

9 —2 9 9 

9eS 


9tt 


U=Uh 


z(u h ) 

Z(u s oo) 


Thus we find 


Sft(cr) = 0 H 


Z (u h ) 


Z(u —> oo) J ’ 

where an is the conductivity evaluated at the horizon and its expression is given by, 


(5.47) 


0 H = 


1 9 XX 


d -1 
2 


g e g g xx 


u=u h 


(5.48) 


where we have used isotropy along the spatial directions (besides u) in the lower dimensional 
theory. Using eq.(5.47), eq.(5.46), eq.(5.43) we get 


»(<t) 


1 

2K 2 <7eff 


d— 1 
9xx 


9x X 


u=u h 


1 

2 s 2 


e * 9zz 


d— 1 
9xx 

9xx 


u=u h 




9xx 


9 zz 
9xx 


( Z(u h ) \ 
y Z(u s oo) y 

/ Z(u h ) \ 
y Z(u s oo) J 

( Z(u h ) V 
y Z(tt —>• oo) y 


(5.49) 


We can now repeat the analysis done in the previous section to evaluate the ratio z^u^oo) k 
by using general coordinate invariance in the underlying higher dimensional theory and 
noting that the gauge field is an off-diagonal component of the metric, eq.(5.39) (for details 

1 With our choice, eq.(5.39), the dimensional reduction results in an action which is not in Einstein 

frame. We could have performed a conformal transformation to bring the lower dimensional action back to 
the Einstein frame. Our end result however will be independent of this choice. 
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see eq.(5.28)). 
Thus we get 


1 


K(a) = — 2 >/G g x £ ^ Z{u h f 
9 XX 


1 


- J—Jr ^ ^ 

- 2K^ WLr 9xx a a 2 

9 XX 9 zz 


u=u h 


u=u h 


u=u h 


The higher dimensional entropy density is 

2tt 


s = ^Vg 


Qxx 


Hence we arrive at the result 


a _ j^D—d 1 Qxx 


u=u h 


(5.50) 


(5.51) 


(5.52) 


S 47T g zz 

Finally, the arguments given in subsection 5.2 allows us to connect r] xz computed in 
the higher dimension to a in the following way 


a 


Thus we find 


9xz — 


Vxz L D ~ d 

S S 


L D ~d 

1 Qxx 
47r g zz 


u=u h 


(5.53) 

(5.54) 


which agrees with the examples we have studied in the previous sections. 


6 Conclusion 

In this paper we have considered a variety of anisotropic examples, and have shown that 
suitable components of the viscosity can become very small in the highly anisotropic case 
and can parametrically violate the bound, eq.(l.l). All our examples have the feature that 
the breaking of rotational invariance is due to an externally imposed forcing function which 
is translationally invariant. E.g. due to linearly varying scalars which give rise to a constant 
forcing function, or due to a spatially constant magnetic held, which was studied earlier in 
[28]. Another common feature in all our examples is that some residual Lorentz symmetry 
survives at zero temperature. In the second half of the paper we show in considerable 
generality that for all cases with these two features, the components of the viscosity tensor, 
which correspond to metric perturbations which carry spin 1 with respect to the unbroken 
Lorentz symmetry, satisfy the relation eq.(1.4). In the anisotropic case the ratio of the 
metric components on the RHS of eq.(1.4) can become very small as T —> 0, resulting in a 
parametrically large violation of the KSS bound. This is indeed true for the examples we 
consider, all of which satisfy eq.(1.4) . 

Besides allowing for a computation of the viscosity with relative ease, the gravitational 
description also provides an intuitive understanding of why such violation of the KSS bound 
may arise. In the absence of isotropy the different metric perturbations break up into 
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components with different values of spin with respect to the remaining Lorentz symmetry. 
Spin 2 components, if present, give rise to viscosity coefficients which satisfy the KSS bound. 
But spin 1 components can violate it. In fact the spin 1 components are akin to gauge fields, 
and the corresponding calculations for these components of the viscosity therefore becomes 
similar to those for conductivity. These are well known in several AdS/CFT examples, and 
also in nature, to sometimes become very small. 

In weakly coupled theories, with well defined quasi particles, we would expect, [37], 
[38], that 

V Imfp 

— n*j - 

s A dB ’ 

where l m fp^dB refer to the mean free path and the de Broglie wave length for the quasi 
particles. This leads to the intuitive expectation that at strong coupling the ratio r]/s ~ 
0(1). However, here we see that at strong coupling, where the gravity description is valid, 
some components of the viscosity tensor in the anisotropic case violate this relation and 
can become parametrically smaller. 

The generality of our result suggests the possibility that this behaviour might happen 
in nature too. It would be very exciting if this can be probed in experiments, perhaps on 
cold atom systems, or in QCD. 

Ordinarily, QCD at finite temperature is described by a homogeneous and isotropic 
phase for which the calculations discussed here are not relevant. This is true even when we 
consider situations which come about due to anisotropic initial conditions, as might arise in 
heavy ion collisions. The behaviour of the QCD fluid in these situations is still governed by 
rotationally invariant Navier Stokes equations with appropriate viscosity coefficients. How¬ 
ever, this could change if a sufficiently big magnetic field is turned on breaking rotational 
invariance 2 . The resulting equilibrium phase could then be highly anisotropic and our 
results, and earlier work, [28], hint that suitable components of the viscosity might become 
small. It has been suggested that such an intense magnetic field might perhaps arise in the 
interior of some highly magnetised neutron stars 3 , see [39], [40] and [41], It has also been 
suggested that strong magnetic fields might actually arise in the highly relativistic heavy 
ion collisions (see [42], [43] and [44]), although in this case the transitory nature of these 
fields must also then be taken into account. 

Turning to cold atom systems, the unitary Fermi gas has also been observed to have 
a value of r\/s close to the KSS bound. Perhaps some way to introduce the breaking of 
rotational invariance can be found in this system. It would then be very interesting to 
examine the resulting behaviour of the viscosity tensor. Even at small anisotropy one 
might be able to see a trend where some components start getting smaller than the bound. 
A natural way to incorporate anisotropy in this case might be to consider the effects of an 
asymmetric trap 4 . 

2 A magnetic field of order 10 16 Tesla or so is needed in order to contribute an energy density comparable 
to the QCD scale ~ 200 Mev. 

3 We thank Gergely Endrodi and Gunnar Bali for a discussion on this issue. 

4 We thank Mohit Randeria for very helpful discussions in this regard and also for his comments about 
the spin diffusion experiments. 
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It is worth mentioning that the spin 1 viscosity components, which become very small 
in onr work, govern the diffusion of the momentum components oriented transverse to 
the direction in which the initial inhomogeneity is set up. For example, take a case with 
anisotropy in the z direction. If the momentum along the x direction, p x , is now taken to 
have an initial gradient along the 2 direction, then its diffusion is governed by the viscosity 
component rj xz , with diffusion length 


D ± 


Vxz 

Hf' 


( 6 . 2 ) 


where s is the entropy density. A small value of — then gives rise to a small value for the 
diffusion constant r ^D± in units of temperature. 

It is perhaps worth mentioning in this context that there have been some recent mea¬ 
surements of spin diffusion in the unitary fermi gas system 5 6 . In three space dimensions, 
with rotational invariance intact, the transverse spin diffusion constant is measured to be 
close to the bound which arises from standard Boltzmann transport theory based on quasi 
particles, see [45]. However, in a quasi-two space dimensions [46], it was found that the 
transverse spin diffusion constant is about three orders of magnitude smaller than this 
bound. It would be worth exploring if these observations can be related to the results 
presented here. 

We have not analysed the stability of the anisotropic solutions discussed in this paper 
in any detail. For the one dilaton case this question was analysed at considerable length 
in [27] and no instabilities were found. This suggests that some examples studied here, 
e.g., the two dilation case, also could be stable. We leave a more detailed analysis of this 
question for the future. It is worth noticing that if an instability appears, it will be when 
the temperature T ~ p, where p is the scale of the anisotropy. As a result one expects 
0( 1) violations of the bound for such systems as well, although not violations where the 
viscosity becomes parametrically small. On a more theoretical note, it would be worth 
obtaining string theory embeddings of the anisotropic systems we have studied here and 
examining if they are stable. Some embeddings for the axion dilaton system were studied in 
[15] and for the one dilaton case in [27] and were found to be unstable, since they contained 
fields which lay below the BF bound of the near horizon geometry. In another instance, 
e.g. [23], though, a stable supersymmetric system with anisotropy was found where suitable 
components of the viscosity become vanishingly small at low temperatures, just as in our 
analysis here. 

We have discussed situations where the breaking of rotational invariance is explicit, 
due to an externally applied source. It would also be interesting to extend this analysis to 
cases where the breaking is spontaneous. Examples are known on the gravity side of such 
phases in the literature, see, e.g., [47-56]. Another direction is to consider Bianchi spaces 
which have been discussed in [57, 58], and which describe homogeneous but anisotropic 

5 The anisotropy force in this case would act in the a direction. This force does not directly enter in the 
diffusion equation for p x . For significant anisotropy, p/T 1, the force is big, and as a result the fluid 
cannot move in the z direction at all. This follows from the bulk geometry, e.g. AdSi x R in the case 
considered in section 2, where Lorentz invariance along the s direction is manifestly broken. 

G We thank Sean Hartnoll for bringing these experiments to our notice. 
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phases in general. Some discussion of transport coefficients in such phases using the gravity 
description can be found in [59]. 
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A Numerical interpolation from the near horizon AdS 3 x R x R to asymp¬ 
totic AdS$, 


Our action consists of gravity, a massless dilaton (j) and a cosmological constant A, in 5 
space time dimensions, 


Sbulk 


1 

2 ^ 


d b xy/^g 


(r + 12A 



(A.l) 


Here 2 k 2 = 16 irG is the gravitational coupling ( G is the Newton’s Constant in 5 dimensions) 
and we set A=l. 

It is easy to show that this system admits an AdS§ solution with metric given by 


ds 2 


— u 2 dt 2 + 



+ u 2 {dx 2 + dy 2 + dz 2 ) , 


(A.2) 


and the dilaton is kept constant. 

We now show that starting with the near horizon geometry given by eq.(4.3), one 
can add a suitable perturbation which grows in the UV such that the solution matches 
asymptotically to AdS§ metric as provided in eq.(A.2). 
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Figure 2. Log-log plot showing the numerical interpolation of near horizon AdS 3 x R x R to 
asymptotic AdS$ , with pi = 1 , p 2 = 1. 


This perturbation is given as follows- 

gttiu ) = 2 u 2 (1 + 5A(u )), 

9uu= 2 u 2 (1 + 5A(u)) 
g X x(u ) = 2u 2 (1 + SA(u )), 

gyyH = y (1 + hC'(u)) 
gzz(u ) = y (1 + 5Z)(w)) 


with 


(A.3) 


5A(u) = ai 6C(u) = C 2 u u , SD(u ) = C 3 u u (A.4) 

ai = -(—5 + 2\/5 )(c 2 + c 3 ), iz = \/5 — 1 . (A. 5) 

5 

The numerical analysis is carried out using NDSolve in mathematica. For the case pi = 1 
, pi = 1 the suitably chosen values for c 2 and c 3 are as follows 


c 2 = 85, c 3 = 85. 


(A. 6 ) 


By adjusting the coefficients c 2 , c 3 to the above values one can ensure that the asymp¬ 
totic behaviour of the metric eq.(A.3) agrees with eq.(A. 2 ) at large u, say u=100000 ; 

The plots in Fig (2) show the metric components as a function of u. These plots were 
obtained by numerical interpolation for the case pi = 1 , p 2 = 1 and c 2 = 85, c 3 = 85 . 


B Ratio of normalizable over non-normalizable mode near boundary 

Here we check that asymptotically the canonical momentum n goes to a constant indepen¬ 
dent of u . To see this , we consider the action 

Sbuik = 2^2 / (r + 12A) . (B.l) 


Numerics 

Fit 
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we get the following solution for AdS$ (setting A=l). 

ds 2 = u 2 dt 2 + + u 2 dx 2 + u 2 dy 2 + u 2 dz 2 ^j . (B.2) 

The metric perturbations go like u 2 ( 1 + ^i) where C\ is constant. 

Hence, using eq.(5.19) and eq.(5.17) we find that 



(B.3) 

Plugging in the higher dimensional metric components from (B.2)we get n(u) = ^Ci 
which is independent of u. Thus asymptotically, the ratio of the normalizable to the non - 
normalizable mode behaves like \C\. 


C Conductivity formula in terms of horizon quantities 


In this appendix , we show the derivation of (5.23) following [35]. The electrical conductivity 
is defined in (5.20) as 

n(tt,o;) , 

The real part can be written as 


(C.l) 


&(<T) = & T 


n(u, a;) 




n(u, uj)Z(u , — uj) 


iuZ(u,Uj) I l«->oo,o ;^0 \iujZ(u,Uj)Z(u,-Uj) J'u-^oo,oj^O 


n(n, uj)Z{u, —uj) 
cuZ 2 (u) 


I U->-CO,UJ—>-0 


O ' ^ n ( n , uj)Z(u, —uj) 
ujZ 2 (u) 


I U —>-00 ,Cd—^0 * 


(C.2) 


Here we used the fact that Z(u,uj) ~ Z(u) is real to leading order when u —> 0. 
We now proceed to show that" 


— ^[n(?z, uj)Z(u, — cn)] = 0, 
du 


(C.3) 


This can be seen as follows 


d 


du 


d 


du 


d 


—9? I N(u) — Z(u, uj) Z (u, —uj) ) = 3 s — N(u) — Z(u,uj) Z{u , —uj) 


-du 


d 


du 


+N (u) — Z (u, uj) — Z (u, —uj) 
du du 


Using (5.16), r.h.s of above equation reduces to 


— M(u)Z(u,uj)Z(u, —uj) + N(u) — Z(u,ui) — Z(u, —uj) 

du du 


(C.4) 


(C.5) 


IT(k,u;) = 5Z ,(S = —7^2N(u)^Z(u,uj) , hence $y[n(u, u>)Z(u, — w)] behaves like a current 
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which is equal to zero since the quantity in the bracket is real. Thus 9[n(w, oj)Z(u — cj)] 
can be evaluated at the horizon i.e. at u = Uh- 

Demanding regularity at the future horizon , we can approximate the behaviour of 
Z(u,u) as follows 

Z ~ e -Mt+r*) 


where r* is the tortoise coordinate, 


r* = 


9uu , 
- au. 

9tt 


Using eq.(5.19) and 


lim — Z(u,u) = — ioj lim li^Hzfu) + 0 (uj 2 ). 

u=tu h dlL u=ru h y gtt 


we get (in the limit oj —> 0) 


(C.6) 


(C.7) 


(C.8) 


\ / U=Uh ' ' 


= cr H 


= A ( x/=s4-s“ 9““ ) ( 7 . z( " ,i) , 

2 k 2 VV 9 « si 


Z(tl|.) V 
Z(u —> oo) i 


(C.9) 


where an is the conductivity evaluated at the horizon and its expression is given by. 

1 


a H = 


2 k 2 g- 


eff 


u=u h 


(C.10) 


where we used isotropy along the spatial directions in the lower dimensional theory. 
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